JATC SC..OOL 
C1\LIF0 R ‘ 1IA 93943 



\ 



NAVAL 



POSTGRADUATE 

Monterey, California 





THESIS 




T2 1 ^ > 3 7 



SECURITY CLASSIFICATION OF THIS PAGE (Whan Dala Sntorod) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4. TITLE Subtitio) 

Computer Programs for Helicopter High Speed 
Flight Analysis 


5. TYPE DF REPORT 4 PERIOD COVERED 

Master's Thesis 
September 1983 


6. PERFORMING ORG. REPORT NUMBER 


7. AUTHOR^ 

Waldo Francisco Carmona 


6. CONTRACT OR GRANT NUMBERf*J 


S. PERFORMING ORGANIZATION NAME ANO AO OR ESS 

Naval Postgraduate School 
Monterey, California 93943 


10. PROGRAM ELEMENT, PROJECT, TASK 
AREA A WORK UNIT NUMBERS 


11. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93943 


12. report DATE 

September 1983 


13. NUMBER OF PAGES 

137 


14. MONITORING AGENCY NAME 9 AODRESSf/f difloroni from Coniroiiing Offico) 


15. SECURITY CLASS, (of thta report) 


15a. DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 



U. DISTRIBUTION STATEMENf (ot thi* Report) 



Approved for Public Release, Distribution Unlimited. 



17. DISTRIBUTION STATEMENT (ot tho obotroct ontorod in Block 20, It difforont from Roport) 



IS. SUPPLEMENTARY NOTES 



IS. KEY WORDS (Conttnuo on row oido it nocoooory ond idontify by biock numbor ) 

Blade Stall 

Compressibility 

Helicopter Computer Program 



20. ABSTRACT (Continuo an rmvoroo §ido it nocoooory ond idondty by biock numbor ) 

This report gives the user of the HP41-CV handheld programmable calculator 
or the IBM 3033 computer, a blade element method for calculating the total 
power required in forward, straight and level high speed flight for an 
isolated rotor. The computer programs consist of a main program which 
calculates the necessary dynamic parameters of the main rotor and several 



00 



FORM 
1 JAN 73 



1473 



EDITION OF 1 NOV 69 IS OBSOLETE 

S/N 0102- LF-014- 6601 



1 



SECURITY CLASSIFICATION OF THIS PAGE (Whon Dmtm Kntormc 



SECURITY CLASSIFICATION OF THIS PAGE (Whm% Dmtm Enfft Q 



subroutines which calculate power required as well as maximum forward 
velocity, stall onset velocity, and velocity for best endurance. 



j 

i 

S N 0102- LF* 014* 6601 



2 SECURITY CLASSIFICATION OF THIS PAGEfP?i«n Dmf Zntmrmd) 



Approved for public release, distribution unlimited 

Computer Programs for Helicopter High Spsei Flight Ar.aly 

by 

Waldo F. Carmona 

n 

Captain, United States Army 
3. A., University of Dayton, 1973 



Submitted in partial fulfillment, of the 
requirements for the degree of 



MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 



from the 



NAVAL POSTGRADUATE SCHOOL 
Sectember 1933 



-V I 








ABSTRACT 

This report gives the user of the HP41-CV har. field 
programmable calculator or the I3M 3033 computer, a blade 
element method for calculating the total power required in 
forward, straight and level high speed flight for an 
isolated rotor. The computer programs consist of a main 
program which calculates the necessary dynamic parameters of 
the main rotor and several subroutines wnich calculate oower 



required as well as maximum forward velccity, stall onset 
velocity, and velccity for best endurance. 
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INTROD UCTION 



A. GENERAL 



The basis for helicopter rotor analysis s as developed in 
the early 1920's when Glauert extended propeller theory to 
the special case of rotating wings. Since that time, the 
development of the digital computer has permitted many 
improvements cn Glauert' s analysis. 

In order to develop a method for predicting the total 
rotor power required in forward flight, it is necessary to 
develop a method that accurately predicts rotor dynamics. 



The prediction of 


rotor 


dynamics in 


forward flight is a 


complex one. Typic 


ally , 


a helicopter 


rotor blade encounters 


a flow environment 


which 


changes rapi 


dly as it moves around 



in azimuth. 

In forward flight, rotor blade sections are subjected to 
azimuthal variations in net only angle of attack but also 
Nach number. As a result, comprehensive performance 
analysis of a helicopter is much more involved than that of 



a conventional aircraft. 

Recent helicopter design trends have been in 
direction of increasing the maximum forward velo 
possible as well as higher blade tip speeds. Sin 



the 



ity 
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speed of the helicopter is added to the speed of rotation of 
the advancing blade, the highest relative velocities occur 
at the tip or the advancing blade. Whan the section Mach 
number of the blade tip exceeds the critical Mach number of 
the airfoil, compressibility effects result. These effects 
include a large increase in profile drag and change in 
pitching moments, therefore creating additional power 
re guirements. 

In all current helicopters there is a tendency for the 
retreating blade to stall. Just as the stall of an airplane 
wing limits the low speed perfomance of an airplane, the 
stall of a rotor limits the high speed potential of a 
helicopter. The relative velocity of the retreating blade 
decreases as forward speed increases. However, the 
retreating blades must produce the same amount of lift as 
the advancing blade. Therefore, as the relative velocity of 
the the retreating blade decreases with forward speed, the 
blade angle of attach must be increased to equalize lift 
throughout the the rotor disc area. As this angle increases 
the blade will eventually stall at some forward speed. For 
these reasons, stall and compressibility thus combine to 

fo a helicopter at higher 

1 1 



limit the performance 



s o ee d s 



For efficient design of a helicopter, the helicopter 
designer should have the analytical tools necessary to 
predict the performance of a helicopter. This report gives 
the user of the HP41-CV and the IBM 3033 computer a means of 
calculating the total power required in forward, straight 
and level flight. 

B. OBJECTIVE 

The objective was to provide an interactive computer 
program which can oe used during the initial design stage to 
estimate the total power required by the main rotor of a 
helicopter, as well as maximum forward velocity possible, 
stall onset velocity, and best endurance velocity. 

Additionally, the desired accuracy chosen at the 
beginning of this effort was to obtain angle of attack 



within 


a 


plus or minus 


one 


half 


degree, 


power e 


srina 


res 


within 




en 


per cent of 


the 


actual 


power 


required, 


and 


to 


obtain 




he 


needed accu 


racy 


in as 


short a 


running 


rime 


as 



oossible. 
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II 



£E2BL2ii 2JSI2I2I2S 

A. DESCRIPTION OF PROBLEM 

Fez a helicopter in level flight, the maximum forward 
speed is limited by the power available as well as stall and 
compressibility effects. It is therefore advantageous to 
develop a simple-to-use computer program that estimates 
blade stall and compressibility effects on helicopter 
performance. 

3. METHOD CF SOLUTION 

In forward flight, the aerodynamre environment of the 
rotor blade varies as the rotor blade rotates with respect 
to the direction of flight. The method chosen to obtain the 
objective accuracy ignores any variable not immediately 
impacting on the accuracy of the determination of angle of 
attack and power requirements. 

The method utilized herein uses a combination of 
momentum and blade element theory to perform rotor 
performance calculations. This theory is initially used to 
determine the induced, profile, and parasite power required. 

The solution method then utilizes a blade element analysis 
to predict the cyclic, collective, and inflow angles 
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3.3 socia red 
the calcui 
az imuthal 
and snail 
attack and 



with the rotor. These angles a 
anion of nhe rotor blade's angle 
positions of 90 and 27 0 degrees, 
power are then estimated as a fan: 
forward velocity. 



r 3 * h 3 n u 5 3 'i 
of arrack an 
Compr esstbil 
ticn of angle 
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DESCRIPTION OF THE PROBLEM SOLUTION 



A. GENERAL 

The calculat icn/prediction of helicopter perforaar.ee is 
primarily a matter of determining the power required and 
power available over ohe desired flight regime. The power 
information may then be used to predice the operational 
capabilities of the aircraft. 

For the ideal helicopter (no losses* , in forward, 
straight and level flight, power required can be subdivided 
into five parts: induced power, P , required to produce 

rotor thrust; profile power, ? 0 , required to overcome the 
shin friction and pressure drag of the main rotor; parasite 
power, ?p, required to overcome the parasite drag of the 
helicopter; compressibility power, P , required to overcome 
tne increase in profile drag when the tip Mach number 
exceeds the drag divergence Mach number of the airfoil; and 
stall power, ? s , required to account for the increase in 
rotor torque as a result of retreating blade stall. The 
HPhl-cv and IBM 3033 computer programs contained herein 
provides a simple, quick means of evaluating the power 
required by an isolated rotor. 
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ASSUMPTIONS 



The major task ip. helicopter performance analysis is the 
determination cf rotor forces, angles, and power. The 
method chosen ignores any variable that does no* directly 
impact on the desired accuracy. Therefore, in order to 
simplify the computational process used, the following 
assumptions are made: 

1. Steady flow through the rotor system. 

2. Small angle approximations are a valid 
representation of the real world phenomena. 

3. All blades considered are rectangular (non- ta perea) 
with only uniform twist possible. 

4. Hinge off-set is zero (i.e., the thrust vector 
passes through the C . G. ) . 

5. The stall angle on most helicopter blades can be 
approximated by the angle that occurs at CLMAX for a 2-D 
airfoil. 

6. The helicopter is trimmed. This implies that the 
sum of ail the moments about the center of gravity (C.G.) is 
zero, ail forces are in balance, and that no lateral 
flapping is present. 
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7. At stall onset, the value of section Iraq 
coefficient jumps approximately 0.08. 

8. The rotor rotates counterclockwise. 

C. NOMENCLATURE 

"Standard" nomenclature is used. Appendix A contains a 
alphabetical list defining all the symbols and mnemonics 
used in the development of this report. Appendix B contain 
an alphabetical list of all HP4 1-CV displays used in the 
HP 4 1 -CV program. 

D. INITIALIZATION 

It is assumed that the initial design of the heliccpte 
has been completed and that the helicopters weight as well 
as the chord, radius, tip velocity, twist, zero- lift drag 
coefficient, and number of blades of ras rotor are known. 
Finally, it is assumed that an initial estimation of 
equivalent flat plate area is available, and that the 
forward velocity of the helicopter is known. 

E. ROTOR DYNAMICS 

There are four rotor parameters whim will help expedit 
later calculations. The first of these is the rotor advanc 
ratio, y ; u is the ratio of the helicopter's forward 
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velocity :o -he rotational velocity. The advance retro can 
be rsoreser.ted as 



(D 



= V v Co s ct 
■ v m 



(dimensionless) 



Applying small angle approximations the advance ratio 
beco mes 



( 2 ) 



y = v, 



T 

where , 



* ( fl R) 



The second dimensionless ratio tha 
calculated is the inflow ratio, X . The 
ratio cf the net velocity up through the 
tip speed. Per the near hover case, y 
and the inflow ratio can be approximated 



t needs to 
inflow ratio 
rotor system 
< 0.1, the T 
by 



be 

is the 
to the 
n a = 0 



(3) A = - /CT/ 2 

In forward flight, the calculation of X requires the 
determination cf the angle of attack of the rotor disc. 
Letting the angle of attack between -he disc plane and the 
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in co ming 


f ree-strea m 


velocity be , 


and assuming angles - o 


be small. 


th.e inflow 


velocity can be 


calculated using 


(3a) 


II 

r< 


C T 


+ u Tan a 3 



2 AW 



where 

a 2 = -Tan (D p /L ) = -Tar. (D p /S) 

D p = (P p * 550)/ V F 

The last dimensionless parameter that needs to be 
calculated is solidity. Solidity is the fraction of the 
rotor disc area than is composed of blades. For a blade cf 
constant chord (i.e. , non-tapered) solidity can be expressed 
as 



(4) a = (b * C)/(?I * R) 

Finally, the tip loss factor, 3 TL , must be considered. 
The tip less factor is used to account for the loss of lift 
that a rotor blade ex parlances due to flow from the rotor's 
lower surface to its upper surface. The tip less factor of 
a rotor can be approximated by 



( 5 ) 



3 = 1.0 - /2 * CT / b 

TXj 
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?. VELOCITY CALCULATIONS 

There are four velocities that are of interest. The 
first of these is the downvash velocity , w. Assuminq 
steady flow through the rotor, the downwash velocity can be 
approximated by 



(6) w 



W 



2 P 




(f“/s) 



NOTE: This equation is not valid for small 

values of forward velocity. 



The second velocity that needs to be calculated is the 
stall onset velocity, V s . The stall onset velocity is the 
velocity at which the retreating blade tip first exceeds the 
static stall angle. The forward spaed at which stall onset 
is first noted can be approximated by the velocity for best 
range, V3R . This is due to the marked increase in profile 
power required at speeds higher than velocity for best 
range. A typical set of power curves for a helicopter are 
shown in Figure 3.1. 

The forward velocity for minimum ?/'/ (i. e. , best range 

velocity) is easily found graphically on the power required 
curve as the point where a straight line through the origin 
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is tangent 
carves are 
needs to be 
best range, 
ex or ession . 



o the curve (see Figure 3.2). Since the power 
not initally known, an analytical expression 
developed which will estimate the velocity for 
The following is the derivation of cns such 
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As shown below 



the parasite power required 



equal no the profile power required, ana one equality is then 
solved for an equation in which forward velocity is the 
variable, the result is a cubic equation ir. which the 
largest rcct defines the point where tne prcfnle power and 
parasite power are equal. 




D 



*0 



1/2PfV p 3= 1/8 a CdopA D (1 + 4.25 y 2 ) 



(7) y- (4.25 /uf) cdo a d v t y - (a/4 f ) cdo x D n 



Next, the largest root of equation 7 needs 



be 



determined. This can be simply done using the cubic root 



equation. Letting 



p = * (4 .25 c/4f) Cdo A d V t 



- ( a /4 f ) Cdo A d V t 3 



Equation 7 can new be written as 



(7a) 




0 



Substituting 
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( X - p/3 ) 



for V F in equation 7a yields an equation of the form 

(7b) X3 + aX + b = 0 

where , 

a = -1/3 ?2 
b = 1/27 (2p2 ♦ 27r) 

The largest root of equation 7b is then given by 

(8) V = A + 3 

w here, 

A 
3 

Equation 8 can now be used as the initial approximation for 
stall onset velocity. 

The third velocity which is of interest is the 
helicopter's maximum forward velocity, VhAX. As shown in 
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Figure 3.2, 
intersection 



the maximum forward velocity is aivan oy -.he 
of the power required and power available 
curves for a given gross weight and altitude [3ef. 1]. 



POWER REQUIRED VS VELOCITY 




FORWARD VELOCITY (KT) 



Figure 3.2. VMAX,V3E, and 7S from Poser Available Curves. 



Therfore, whenever V? > VXAX , there is 
available to sustain level flight. 

The power-limited maximum speed may 



i nsuf f 



b e est 



cient oower 



mated by 



24 



(9) 



MAX 



1/3 

“of" ^ AVAIL* ? I ' P o } (ft/ £9 



Equation 9 can be simplified by assuming chat the 
required at maximum speed is about the same as che- 
at hover. Therefore, assuming 



( 10 ) 
and that 



(? ~? - ?) = (? - D ! 
' AVAIL I o' ' HOVER ‘ o ' 



( F - ? ) = ( P ) 

AVAIL o I HOVER . / 2 P A 



V" 



equation 9 can then be written as 



(10a) V iMAX = * (Vf/A D ) 



1/3 



(ft/s) 



where , 



(10b) 






VH 



P K 



The last velocity that needs to be calculated 
best endurance velocity, VBE. In the normal operati 
the total helicopter power can be represented by 



( 11 ) 



? T = ? I * P o + P P 



C) 

?c we r 
r e q u ir ed 



is the 
ng ranee 
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Assuming than the variation of profile power with forward 
velocity is negligible, the velocity for best endurance 
(also the best rata of climb velocity) ran be found. If 
eguaticn is differentiated with respect to forward 
velocity, V F, and' is set equal to zero, it can be seen that 



(11 a) P z = 3 ? p 

or 

(11b) W2 = 3 p f V p 2 

2 o Vbe “ 



Solving equation for the best endurance velocity, 733, 
results in 



(11c) 




(f-/s) 



G. INITIAL POWER CALCULATIONS 

As forward velocity increases, the induced power 
decreases, the profile power increases slightly, and the 
parasite power increases until it becomes the dominant loss 
at high speeds [Ref. 2]. For forward, straight and level 
flight, the induced power can be calculated by 
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( 12 ) 




If tip losses are taken into effect (equation 5), the 
induced power now becomes 



(12a) ? x = ( 1/3^0 * 

(TL) 

Since the induced power required in ground effect is less 
than that required out of ground effect, equation 12 should 
be written as 



(12b) P T = (V^) * (G2) * ?! 

(TL+GE) 
where , 

GE = (-0. 1276 * (h/D) 4 + 0.703 * (h/D) 3 

-1 .4569 * (h/D) 2 + 1.3432 * (h/D) 

+ 0.5147) 



The profile power required is given by 



(13) P Q = acdop (1 + 4.25 y 2 5 

4400 



(hp) 



Finally, parasite power given by 
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( 14 ) 



?P = Pf Y 
1100 



(hp) 



H. ANGLE CF ATTACK 

The determination of angle of attar* at the azimuthal 



positions cf 90 


and 270 degree 


s are important in 


the 


determination cf 


compressibility 


and stall effects. 


T - 


because cf this 


effect that the 


dynamics of the blade 


310 * 


are important in 


analyzing helic 


opter performance. 





The angle of attack of the rotor is a function of 



radius, r, and azimuthal position, \p 
illustrates the sign convention used in 
angle of attack. The angle of attack cf 
estimated by, 



. Figur e 3 . 3 
deter min in g blade 
a rotor can thus be 



(15) a (r , W ) = 0 - VpBcos ct + w - 7 a 3 

QR + VpSin ip (rads) 

w h er e , 

(15a) 0 = 0 Q + 0 T + 0Cos ? + 0Sir.Y + K B 6 

(15b) B = B q - a 1 Cos ip - bl Simj) 
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ANGLE OF ATTACK 




Figure 3.3. Rotor Angles in Longitudinal Plans. 

To determine angls of attack, in is then necessary to 
determine the longitudinal collective and cyclic angles 0 
and © 2 . This can be accomplished by expressing the 
coefficient of thrust as 
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(16) 



C T = _|£ ATI + (0 o + K b 8)T2 f0 T T3 + (0 2 - K b d1)T 

wher e , 

T 1 = .5(32 + , 5 P 2 ) T3 = . 25B 2 < B 2 ♦ y 2 > 

TL TL TL 

T2 = (. 523 + . 5 u 2 3 ) 74 = . 5 V (3 2 + .25y 2 > 

TL TL TL 

Additionally, the longitudinal flapping coefficient al, 
needs to be determined. The longitudinal flapping 
coefficient can be written as 



(17) al = A Al 1 + ( 0 + K_e ) A12 + 0 A13 + (0 - K b 1 ) A14 

o^o T 2 B 



where , 

D 1 = (3 2 - . 5 U 2 > 
TL 



All = .25 (.5 ^ 2 L - y3 /3 ) A 13 = 2 y B^ 2 

3 2 3 1 Dl 

TL 



Al 2 = 8 y H__ 

TL 

3 Dl 



A 1 4 = 5 2 +■ 1 . 5 y 2 

TL 
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Assuming that there is nc lateral flapping ana that the 
effect is zero (i.e., Kg = 0.3), equations 16 and 17 beccm 
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( 18 ) 



2 C T = ATI +0 T 2 + 0 T3 +0 T4 

a o T 2 



(19) al = A All + 0 A 1 2 + 0™,A 1 3 +0 A14 

0^2 



NOTE: Since the analysis is only locking an 

the azimuthal positions or 90 and 270 
degrees, in car. be seen from eguation 
156 that the contribution of tne 
longitudinal flapping coefficient, al, 
wilx always be zero at these positions. 



Equations 18 and 19 new represent a set of simultaneous 
equations in which the only unknowns are the collective and 
cyclic angles and can thus be determined. Knowing the 
values of the cyclic and collective angles, the blade tip 
angle of attack at the 90 and 270 degree positions can be 
estimated by 

(20, a 9Q = e o + e T - e 2 

(21) a 270 = • ® 2 + e T 

NOTE: The anqle of attack, a , is defined 

oositively if the disk olane is nose 
ud, see Figure 3.2. 



+ A 
l+y 



+ _A_ 

1+ U 
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STALL POSSE 



In the helicopter stall nor 
the retreating blade, sines the 
usually at the blade. tip. As t 
the stalled area of the rotor bl 
At the higher values of u, t 
required are great and therefor 
Assuming a jump of 0.08 in the v 
that the rotor area within which 
segment of minimum dimensionlss 
stall area is symmetric about $ 
New found that the effects of ti 
the higher values of u are large 
high speed flight by 

(22) C ps = <1 - p 



maliy starts at the tip of 
highest angles of attack are 
he forward speed increases, 
ade spreads inboard, 
he effects of stall on power 
e need to be estimated, 
alu e of Cdo at stall onset , 
the blade stall exists is a 
s radius X s and that the 
= 270 degrees, Castles and 
p stall on power required at 
and can be approximated for 

) 2 (1 - Xg) /I - X| 



where X s is the nondimens icnal radius outboard of which the 
retreating blade is stalled [Ref. 3], The dimensionless 
radius, X g , can be estimated by equating the section angle 
of attack, at 'i> = 27 0 degrees, to AN AX 'Ref. 4]. Setting 
equation 15 at ij ; = 270 degrees equal to A MAX results in the 
quadratic listed below. 
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(23) 


(X -y) = © q (X -y) t 0 T X 2- © T Xy-0 2 (X -y) + A +x al 


Applying 


the quadratic formula to equation 23 yields the 



roots 



(23a) 


x s = ~ 3 s * l/B s * 40 t c s 

20 t 


(23b) 


X o * - B S - /B I - “ e T C S 
2 0 t 

where , 


(23c) 


T = -0+0 

MAX o 2 


(23b) 


Cg = pT+A 


(23 e) 


5 = -y0 T -r 


Eauat 


:icn 22 is satisfactory for most cases. It is 


do ss ible 


however, for the blade section angle of attack to 


be higher 


: inboard than at the tip creating a situation which 



is ussually refared to as inboard stalling [ Ref . 5], 

For the special case of inboard stalling, see Figure 
3.4, the incremental stall power coefficient defined by 
equation 22 is too large and needs to be corrected. 
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Assuming the stalled region is diamond shaped as shown in 
Figure 3.5 and that the stalled area is symmetric about ip 

3tt/2, it can be seen that as X approaches X , the 

o'' S 

correction to the incremental stall power coefficient, C ps 



3U 



must vanish 



STALL REGIONS 




Figure 3.5. Approximated Stall Area. 

Similiarly, as the average of X and X aoproaches unity, 

o S 

the value cf C ps goes to the value defined by equation 22. 

Therefore, in order to correct for the possibility cf 
inboard stalling, the correction factor, k , is defined such 
that 
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(24) 



k 



1.0 



for - 3 



S 




(25) 



k = B /20 + X 0 
.q t s 



tor -3 



S. < 1.0 




20 



2 



Ths corrected equation for stall power required thus becomes 



J. COM? H ESSI3ILITY POWER 

The individual effects of stall and compressibility on 
rotor profile power are substantial an high aivancs ratios. 
When both effects are present, the losses due to each source 
are difficult to distinguish. Therefore, as a helicopter’s 
forward velocity and tip speed increase, the need for a 
simple estimate of how the compressibility of the air 
influences the rotor performance is neccessary. 

In forward flight, the Mach number of the advancing 
blade is given by 
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(27) 






\l ? (x ♦ y Sir. ip ) 



where , 

( 27a ) » TIP - fiR = ? T 

a sv a sv 

(27b) x = r/R = nondimensional radius 



Since t 


he 


highest 


M a ch 


number occurs 


at the tip of the 


ad va nci 


ng 


blade at 


ip 


90 degrees. 


eg nation 27 can be 



written as 

(28) a * m ( i + y ) 

90 TIP 

Reference 3.5 showed chat the critical Mach number for drag 
divergence can be estimated by 

(29) M = MC3 0 -ml * a * a 

CRIT 90 

Gesssow and Crim, in their investigation of 
compressibility found that the compressibility effect on 
rotor performance was a rapid increase in the profile power 
when the tip Mach number exceeded the critical Mach number 
for drag divergence [Ref. 6]. The increase in profile power 
coefficient due to Mach effects can be estimated by 
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( 30 ) 



( 31 ) 



c pc = 0.0 1 AM d + 0 . 1 ( A^l 3 

w h sr 9 t 



A = K 



90 



1 

GRIT 



0.06 
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IV 



CDNCLUSICNS AND RECOMMENDATIONS 



The objective of ~his project was to provide at easy re 
use, int er-active computer program which could be used 
during the initial design phase to estimate the coral power 
required by the main rotor of a helicopter. The computer 
programs contained herein provide results which are well 
within the objective accuracies (see Figures C.1, C.2, C.3) 

and provide acceptable results as a first cut estimate of 
compressibility and stall power requirements. 

In the development of the computational model, many 
simplifying assumptions were made to ease the amount of 
computation required. The assumption which most impacts on 
the accuracy of the program is that of steady flow through 
the rotor. 

The flow environment encountered by the rotor changes 
rapidly due to the rate of change of blade angle of attack. 
Additionally, rotor operation at high advance ratios also 
produces considerable radial flow along the blade span . 
The steady flow assumption ignores the time-variant aspect 
of rotor aerodynamics. The dynamic nature of the rotor, 
especially when operating at or near the stall regime, 
requires the application of unsteady aerodynamics, and a 
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close examination of 


ho w the 


pi tching 


aomsnts ge: 


the retreating 


blad e 


s t ailing 


af f acts 


controiabilii y [Hef. 


7]. It 


is rher 


efore recoin 



additional investigations consider how the u.is 
aerodynamics, and pitching moments generated, ir.fl 1 
performance of the helicopter. 



a red by 

e d than 
ready 
er.ce th 
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A PPEND ! X A 
NO MEN CLATURE 



Ter a Mnemonic 

a a 

A AD 

ct 3 Al pha3 

a MAX A M A X 

a A9 0 

90 

a A 270 

270 

a 1 al 

All All 

A 1 2 A 1 2 



Def i ninio n 
Slope of airfoil 
section lift curve. 

So tor disc area. 

Disk plane angle of 
attack. 

The steady flow stall 
angle of the airfoil 
(given by a * CL MAX ) . 
Angle of attack of 
the advancing blade 
at if) = 90 degrees . 

Angle of attack cf 
the retreating blade 
it t = 270 deg rees. 
Longitudinal flapping 
Term in definition of 
THETA2. 

Term in definition of 
THETA 2. 



Unit- 



rad 



rad 



rad 



rad 



_ a d 



rad 

dimensionless 



dimensionless 
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A 1 3 A 1 3 

A 1 4 A 1 4 

b b 

S 3STA0 

o 

3 S 3 5 

b tl btl 

c c 

Cdc CDO 

cl max clmax 

Cp C CPC 



C PS 

c s 



CPS 



cs 



Term in defnr.it on cf 
THETA2. 

Term in definition of 
I HETA2 . 

Number of blades on 
the main rotor . 

Main rotcr coning 
angle. 

Term in definition of 
X S * 

Tip loss factor. 



dimensionless 



dimensionless 



dimensionless 



rad 



dimensionless 



dimensionless 



Sean chord of main ft 

rotor blade. 

Main rotor coefficient dimensionless 

of drag at zero lift. 

Maximum coeficient of dimensionless 

lift (2-D). 

Correction to power dimensionless 

coefficient due to 

com pr essibilit y effects. 

Correction to power dimensionless 

coefficient due to stall. 

Term in definition of X g . dimensionless 
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C T CT 

Dp D 

D DP 

Ail D d:id 

6 DEL 3 

3 



T GAMMA 

G. £. G S 
h H 

hp HP 

K KB 



Coefficient of thrust 
of the main rotor. 

Rotor disc diameter. 
Parasite drag of the 
helicopter. 

Term in the definition 
of compressibility 
p o w er . 

Rate of change of blade 
pitch with respect to 
blade flapping. 
Equivalent flat plate 
area of the helicopter 
in forward flight. 

Term in definition of 
3 an d C 

ground effect ratio. 
Height of main rotor 
above the ground. 
Horsepower, 
effect . 



dimensionless 



h 



dime rsi or.less 



dimensionless 



dimensionless 



Lime nsionless 



hp 

dimensionless 
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A LAMB 


Ratio of the net dimensionless 

velocity up through 
the rotor system to 
the tip speed. 


M MACH 


Mach number of rotor dimensionless 

blade. 


m crit mcrit 


Critical Mach number dimensionless 

of advancing blade 
at ip = 90 degrees. 


M m MT 

T 


rip Mach number of dimensionless 

rotor balde. 


M MCRTO 

CRTO 


Critical Mach number dimensionless 

for Cl * 0. 0 . 


V 


Main rotor advance dimensionless 

r a t io . 


ml Ml 


Constant in definition ft-lb 

of critical Mach number 
(ml = 0.113). 


OMEGA 


Rotational velocity. rad/s 


? C ?H 


Power required due to hp 

compressibility effects. 


O D T 

* I 


Induced power. hp 
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D 

* 0 


PO 


Pro fils power. 


hp 


"5 

‘ P 


?? 


Parasite Power. 


h p 


?s 


PS 


Power required due to 


hp 






stall affects. 




H 


R 


Main rotor radius. 


f - 


P 


RHO 


Air density. 


slug/ft 


a 


SD 


Main rotor solidify. 


dimensio n 1 9 s s 


0 


THETA 


Ratio of ambient temp- 


dimensionless 






erature to standard sea 








level temperature. 




0 


THETAO 


Main rotor collective 


rad 


0 




pitch. 








Main rotor lateral 


rad 






pitch. 




0 


THETA2 


Main rotor longitudinal 


rad 


2 










cyclic pitch. 




T1 


T 1 


Term in definition of 


dimensionless 






THETAO . 




T2 


T 2 


Term in definition of 


dimension! ess 






THETAO . 




T3 


T 3 


Term in definition of 


dimensionless 






THETAO . 




14 


T4 


Term in definition of 


dimensionless 
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THETAO 



V I 


71 


Induced velocity. 


ft/s 


7 f 


7 F 


Aircraft forward speed. 


f t/s 


7 


7KT 


Aircraft forward speed. 


kt 


7MAX 


7M A X 


Maximum forward velocity 


kt 



possible. 

7 7 5 Snail onset velocity kr 

u 

(velocity at which A270 
equal A MAX) 



V T VT 

a S7EL 

SV 

v DW 

A A 

x s XS 



X xo 



Rotor tip speed. ft/s 

Speed of sound. ft/s 



Rotor dcwnwash velocity. 
Aircraft gross weight. 
Radius cuxboard of which 
the main rotor is stalled. 
Radius inboard of which 
rotor blade stall may 
be present due no inflow 
ratio and blade twisr. 



ft/s 

lb 

dimensionless 



dime nsio nless 
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APPENDIX B 



22 '±lz£l 2 bog a Ail documentation 

A. GENERAL 

The HP 41-CV program uses 45 program storage registers. 
Prior ro program initialization , the handheld calculator 
should be sized to 46. 

The computer program consists of a main program which 
calculates the necessary dynamic parameters of the main 
rotor and several subroutines which calculate power required 
as well as the maximum forward velocity, stall onset 

velocity, and best endurance velocity. 

The documentation for the HP 41-CV program is device! 
into the following sections: 

1. PURPOSE 

This section describes the intended purpose of 
the program or subroutine. 

2. ASSUMPTIONS 

This section lists any assumptions made which are 
applicable to the program cr subroutine. 

2. EQUATIONS 

This section lists the equations utilized within 
the main program cr subroutine. The primary references for 
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the 



equat ions 



used 



are Aer ofl v n an ics and V/STOL Flight, 



Reference 3.1, Re l ie op ter Theory, Reference 3.2, and 

Ai rcraf t Performance, Reference 3.3. 

4. FLOWCHART 

Both the handheld computer and the I3H 3033 
normally execute instructions in a program in a sequential 
manner unless it is instructed to do otherwise. This 
section will graphically represent the step by step method 
used tc solve the problem as well as the flow cf control 
between the various parts of the program. In a flowchart, 
different types cf operations are indicated by different 
shaped boxes as illustrated below: 




R sctangl e 



For a calculation or 
process other than a 
decision. 
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am 



?du input or output. 



Small Circle 




For an on page 
oo a n action when a 
f 1 o wchart continues 
on the same page, or 
when it is difficult 
to connect two boxes. 




For connection when 
a flowchart continues 
to another page. 
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5. PROGRAMS AND SUBROUTINES 

This section lists all the pro grass or 
subroutines which must be present in the HP ui-cv prior to 
program execution. 

6. PROGRAM LISTING 

This section contains the HP 4 1 -CV listing of the 
program or subroutine. 

B. QUICK REFERENCE TABLES 

The tables in this section will be useful to the user of 
the HP 41-CV handheld calculator as a source of quick 
reference. Table 1 is an alphabetical listing of ail 
calculator displays with an explanation of their respective 
meaning. Table 2 lists all storage registers used and 
describes how they are utilized. 



50 



1 



• 32iflz£V Disol av 5 



TABLE I 



DI 5? LAY 


DEFINIT 10 M 


UNITS 


a = ? 


Slops of airfoil lift curve. 


rad 


b = ? 


Number of main rotor blades. 


N/A 


c = ? 


Chord length of main rotor. 


C „ 


DA = ? 


Density altitude. 


£ • 


Cdc = ? 


Main rotor coefficient of drag 


N/A 




at zero lift. 




FP A = ? 


Equivalent horizontal flat plate 


f “ 




area of the helicopter. 




h = ? 


Height of the main rotor above 


C im 




the ground. 




SC RO =? 


Critical Mach number for the 


N/A 




coefficient of lift equal to 






zero. 




a = ? 


Main rotor radius. 


f j- 


TWIST = ? 


Geometric twist of the rotor. 


rad 


7? (K T) = ? 


Forward velocity. 


Kt 


7T = ? 


Fetor tip speed. 


ft/s 


W = ? 


Aircraft gross weignt. 


lbs 


CD LL = 


Main rotor collective pitch. 
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CYCLIC = 


Main rotor cyclic pinch. rad 


PI = 


Indue ad power required hp 

compensated for tip losses 
and ground effect. 


PM = 


Power required due to hp 

compressibility effects. 


PO = 


Profile power required. hp 


PS = 


Power required due no one hp 

retreating blade stalling. 


V3B = 


Maximum endurance velocity. Kt 


X 

> 

A 

> 


The forward velocity that has M/A 

been input is larger than the 

one calculated by subroutine VMAX. 


7 M A X = 


Maximum forward velocity. Kt 


VS = 


Initial estimate for finding Kt 

onset velocity. 


XS = 


Radius outboard of whicn the M/A 

main rotor is stalled. 


XO = 


Radius inboard of which rotor M/A 

blade stall may be present. 


a9 OD = 


Angle of attack of the advancing deg 

blade at = 90 degrees. 


a2 70 = 


Angle of attack of the retreating dec 
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2 . 



2. H?41 


blade at i|> = 270 degrees. 
-CV Reqistsr Utilization 




TA3L5 II 


STORAGE 


REGISTER 


STORED QUANTITY/USE 


00 


Ground effect ratio (GE) . 


01 


Maximum forward velocity (VMAX) . 


02 


Stall onset velocity (VS) . 


05 


Rotor radius (R) . 


06 


Number of blades ( b) . 


07 


Zerc-lift drag coefficient (Cdo) . 


09 


Rotor height above the ground (h) 


10 


Aircraft gross weight (W) . 


1 1 


Air density (RHO) . 


12 


Lift curve slope (a) . 


13 


Rotor tip velocity (VT) . 


14 


Coefficient of thrust (CT) . 


15 


Tip-loss factor (STL). 


16 


Main rotor induced power (PI) . 


17 


Rotor height to rotor diameter ra 


18 


Compressibility power (PC). 



o (h/d) 
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19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

41 

42 

43 

44 

45 

HE 



Solid! tv (SD). 

Induced velocity (VI). 

Main rotor profile power (PO). 

Advance ratio (HU) . 

Stall power (PS). 

Advance ratio squared (HU 2 ) . 

Forward velocity ( V?) . 

Equivalent horizontal flat plate area (??A) . 
Maximum 2-D lift coefficient (CLMAX) . 

Main rotor parasite power (?P) . 

Main rotor geometric twist (TWIST). 

Angle of attack at 270 degrees (A270) . 

Angle of attack at 90 degrees (A90) . 

Sonic velocity (SVEL). 

Critical Mach number for coefficient of lift 
equal to zero (MCRO). 

Density altitude. 

Scratch pad calculations. 
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C. PROG RA II DOCUMENTATION 

This section contains the necessary documentation, for 
the H? 41-CV computer program. The main program as well as 
all of the subroutines used in the solution of the problem 
are outlined in this section. 
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1 . 



21 ait £ro ar^a 

(a) PURPO S3 

This program calculates the dynamic parameters of 
the main rotor which are necessary for calculating the 
total main rctor power required in forward, high speed 
straight and level flight. It additionally controls the 
execution sequence of the various subroutines which are 
used to calculate main rotcr power required, in terms 
of horsepower, as well as maximum forward velocity, stall 
onset velocity and velocity for best endurance in knots. 

(b) ASSUMPTIONS 

(1) All angles are small. 

(2) Steady flow through the rotor. 

(3) All rotor blades are rectangular (non-ta pered) 
with only uniform twist being possible. 

(4) Only the first hamonic of flapping is necessary 
for calculating power required. 

(4) The effective dimensionless radius can be 
apprcxima t eay by the tip-loss factor. 

(6) The thrust vector passes through the C.G.. 

(7) The static stall angle fcr blades on most 

helicopters can be approximated by the angle at 

which CLMAX occurs for the 2-Q airfoil. 
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( 3 ) 



) 2QU A 



The 6 3 affect, or the result of cocking 
flapping axis of the blade sc : ha: ins pit 
varies as the blade flaps is zero (i.a., n 
lateral flapping is present). 



TIOHS 



u = v p /v T 

V (ft/s) = V * 1.68394 
F 

T 1 = .5 (E^£ + .5 U 2) 

T 2 = . 3 34B 3 ♦ . 5 y 2 B 

TL TL 

T3 = . 25B m 2< 3 2 + y 2 ) 

TL TL 

T 4 = . 5y (9 2 + . 25 y 2 > 

TL 

D1 = (B tl 2 - .5 y 2 > 

All = 4 (.5 B 2 - yJ / 8 ) /D 1 
TL 

A 12 = 8yB TL /3D1 

A 13 = 2yB T] ? / D 1 

A 14 = (B 2 + 1 . 5 y2> / D 1 
TL 

(2CT/aa) = ATI + 0 q T2 ♦ 9 T T3 + 9^4 

0 = AA1 1 3 A 12 + 9_ A 1 3 + 9 A14 

o T 2 

Amax = CLMAX/CLA 



a27 0 = 9 - 9 

o 2 

r = AMAX -0 

o 

C s = yr + A 

b s * - u« T - r 

*s = 'V <V 

\ = 'V < 3 S* 



+ 9 + ( A / 1 + y ) 

-0 

2 



- 49 t C s )-s / 29 t 

- 49 t C s )-s / 29 t 
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a. FLOWCHART 




i 

B 

V 
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1 

0 



\/ 
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0 



\/ 




e. PROGRAMS AND SUBROUTINES USED 



SD " 


n 


CT " 


PIT " 


n 


LAMB 


POl " 


« 


CPS H 


PP1 " 


tt 


CPC " 


VS " 


n 


a90 ” 


VMAX " 
DEN " 


n 


CNG " 



60 



f. PROGRAM LISTING 



01+LSL "MBS" 



82 


*R=?* 


83 


PROMPT 


84 


STO 85 


85 


-b=?’ 


86 


PROMPT 


87 


STO 86 


88 


•O?* 


89 


PROMPT 


18 


STO 84 


11 


"CdO=?- 


12 


PROMPT 


13 


STO 87 


14 


-«=?■ 


15 


PROMPT 


16 


STO 18 


17 


■VT=’' 


18 


PROMPT 


19 


STO 13 


28 -VF<KT>=?' 


21 


PROMPT 


22 


1.68894 


23 


* 


24 


STO 25 


25 


-FPP=?- 


26 


PROMPT 


27 


STO 26 


28 


-H=?- 


29 


PROMPT 


38 


STO 89 


31 


•a=?" 


32 


PROMPT 


33 


STO 12 


34 ' 


"CLMflX=‘ ? 


35 


PROMPT 


36 


STO 27 


37 ' 


■T«IST=? 


38 


PROMPT 


39 


STO 29 


48 


*MCR8=? 


41 


PROMPT 


42 


STO 44 



43 "Dfi=?* 

44 PROMPT 

45 STO 45 

46*LBL *QGN“ 

47 RCL 25 

48 PCI 13 

49 / 

58 STO 22 

51 Xt2 

52 STO 24 

53 XEQ -SB- 
54 XEQ -DEH- 

55 FS? 03 
58 GTO 85 

57 XEQ -PIT- 

58 XEQ -PQ1- 

59*L8L 05 
00 XEQ -PP1- 

61 flBV 

62 PS? 03 

63 GTO 06 
64 XEQ -VBE- 

65 XEQ -V«OX- 
66 QBV 

67*181 96 

68 XEQ -CT- 

69 RCL 15 
78 Xf2 

71 STO 38 

72 EHTER.t 

73 RCL 24 

74 2 

75 / 

76 + 

77 STO 31 

78 2 

79 / 

88 STO 33 
81 RCL 15 



82 3 
33 YtX 

84 3 

85 / 

86 EHTERt 

87 RCL 24 

88 RCL 15 

89 * 

98 2 

91 / 

92 * 

93 STO 34 

94 RCL 38 

95 RCL 24 
% + 

97 RCL 38 

98 * 

99 4 
188 / 

181 STO 35 

182 RCL 38 

183 RCL 24 

184 4 

185 / 

186 ♦ 

187 RCL 22 

188 * 

189 2 
118 

111 STO 36 

112 RCL 22 

113 RCL 38 

114 * 

115 2 

116 / 

117 RCL 22 

118 3 

119 YtX 
128 8 
121 / 

122 - 
123 4 
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124 * 

125 RCL 38 

126 ENTERt 

127 RCL 24 




131 STO 32 

132 RCL 38 

133 * 

134 / 

135 STO 37 

136 RCL 22 

137 RCL 15 

138 * 

139 8 
148 * 

141 3 

142 f 

143 RCL 32 

144 / 

145 STO 38 
148 RCL 22 
H7 2 

148 * 

149 RCL 38 
158 * 

151 RCL 32 

152 / 

153 STO 39 

154 RCL 38 

155 ENTERt 

156 RCL 24 

157 1.5 

158 * 

159 + 

168 RCL 32 

161 / 

162 STO 48 

163 RCL 14 
if, 4 ? 

165 * 



207 - 
268 / 

209 STO 32 
218 FS? 83 

211 GTO 87 

212 P-3 

213 “CYCLIC=" 

214 fiRCL X 

215 fiVIFy 

216 STOP 

217+LBL 87 

218 RCL 31 

219 ENTERt 
229 RCL 36 

221 RCL 32 

222 * 

007 _ 

224 RCL 34 

225 / 

226 STO 33 

227 FS? 83 

228 GTO 88 

229 R-l! 

238 -COLL=* 

231 PRCL X 

232 PVIEH 

233 STOP 

234 QBV 

235*LBL 88 

236 RCL 33 

237 RCL 32 

238 - 

239 RCL 29 
248 + 

24! ENTERt 

242 RCL 83 

243 1 

244 RCL 22 

245 + 



166 RCL 12 

167 / 

168 RCL 19 

169 / 

178 STO 42 

71 XE8 -LfiNB 

172 RCL 83 

173 RCL 33 

174 * 

175 CHS 

176 RCL 42 

177 + 

178 RCL 29 

179 RCL 35 

188 * 

181 - 

182 STO 31 

383 RCL 83 

184 CHS 

185 RCL 37 

186 * 

187 RCL 29 

188 RCL 39 

189 * 

198 - 

191 STO 32 

192 RCL 34 

193 * 

194 ENTERt 

195 RCL 31 

196 RCL 38 

197 * 

198 - 

199 RCL 34 

288 ENTERt 

281 RCL 40 

282 * 

283 RCL 36 

284 ENTERt 

285 RCL 38 

286 * 



52 



246 / 

247 + 

248 STO 41 

249 FS? 83 

250 GTQ 89 

251 P-B 

252 ■i278=‘ 

253 fiRCL X 

254 OVIEU 

255 STOP 

256 fiBV 

257 XEQ -498- 

258 fiBV 

259 PCI 27 

260 RCL 12 

261 / 

262 RCL 33 

263 - 

264 RCL 32 

265 + 

266 STO 38 

267 RCL 22 

268 CHS 

269 RCL 29 
278 * 

271 RCL 38 

272 - 

273 STO 31 

274 RCL 38 

275 RCL 22 

276 * 

277 RCL 83 

278 + 

279 STO 35 
288 4 

281 CHS 

282 * 

283 RCL 29 

284 * 

285 STO 35 

286 RCL 31 



287 Xt2 

288 STO 36 

289 EHTERt 
298 RCL 35 

291 CHS 

292 x>y? 

293 SF 01 

294 FS? 81 

295 GTO 83 

296 RCL 36 

297 RCL 35 

298 + 

299 SORT 
388 STO 38 

301 RCL 31 

302 CHS 

383 + 

384 RCL 29 
305 / 

386 2 

307 / 

388 STO 38 
309 ■XS=' 
318 fiRCL X 

311 flVIEH 

312 STOP 

313 RCL 3! 

314 CHS 

315 RCL 38 

316 - 

317 RCL 29 

318 / 



321 STO 88 

322 'XO=* 

323 fiRCL X 

324 fiVIEH 

325 STOP 

326 fiBV 



327*LBL 83 
328 XEQ "CPS" 
329 fiBV 

338 CF 81 

331 CF 02 

332 CF 83 
333 XEQ 'CPC- 

334 fiBV 

335 fiBV 

336 XEQ -PT1- 

337 XEQ -CHG- 

338*L8L 89 

339 XEQ 'VS' 
348 ENB 
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Solidiiv 



a. PURPOSE 

This subroutine calculates tne ratio of total blade area 
to the total rotor disc area. 

b. ASSUMPTIONS 
None 

c. EQUATIONS 



b * c * R 



b * c 



PI * R 



PI * R 



d. FLOWCHART 
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PROGRAMS AND SUBROUTINES USED 



y . 



’ Y8S ’ 

f . PROGRAM LISTING 



81*LBL "SB” 
92 RCL 86 

83 RCL 8* 

84 * 

85 RCL 85 
@6 f 

87 PI 

88 / 

89 STO 19 
18 EHD 
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3. Down wash Vs loci tv 



a. PURPOSE 

The purpose of this subroutine Is to compute the Induced 
velocity of a rotor system. 

b. ASSUMPTIONS 

Steady flow through the rotor system. 

c. EQUATIONS 

u 

2 ’ A o V r 

d. FLOWCHART 






start 



i 


f 


calculate 

downwash 

velocity 


3 




store DW 
In R32 


1 


! 



return 



e. PROGRAMS AND SUBROUTINES USED 

" WBS" 
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f. PROGRAM LISTING 



01*1 


LSI 


-H" 


02 


RCL 


- c 

CJ 


03 


; s !=8? 




84 


GTO 


85 


85 


RCL 


10 


86 


0 




87 


RCL 


11 


88 


* 




89 


RCL 


25 


18 


* 




4 \ 


PI 




10 

iu 


* 




17 


RCL 


85 


14 


Xt2 




15 


* 




16 


/ 




i *? 

i ?' 


STO 


64 


18 


GTO 


86 



i?*L3L 85 
28 SF 83 

2i*LBL 86 
22 EHB 
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4 



. Coefficient of Thrust 



a. PURPOSE 

This subroutine calculates the coefficient of thrust for a 
arbitrary rotor. 



b. ASSUMPTIONS 

Steady flow through the rotor system. 

c. EQUATIONS 



d. FLOWCHART 



start 



* 


l 


calculate 
coeff of 
thrust 




r 


store In 
R14 


> 


< 



return 
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e. PROGRAMS and subroutines used 



‘ WBS * 

f. PROGRAM LISTING 



01+LBL "CT" 
82 RCL 95 

93 Xf2 

94 PI 

95 * 

96 RCL it 

97 * 

88 RCL 13 
09 Xt2 
18 * 

11 1/X 

12 RCL 19 

13 * 

14 STO 14 

15 END 
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5. Induced Powe 



(a) PURPOSE 



This subroutine calculates the power required by the rotor 



to produce thrust at. hover and forward flight. Additionally, 
this subroutine corrects for tip losses (losses in lift at the 
tips due tip votices) as well as ground effect. 

(b) ASSUMPTIONS 

Steady flow through the rotor system. 

(C) EQUATIONS 



PI 




I 



B 



1.0 - V 2 CT 



D 



GE = -0. 1276 (h/D) * + 0708 (h/D) 3 



- 1. 4569 (h/D) 2+ 1 343 2 (h/D) +0.5147 
PIT = ( 1/3 ) * (GE) *PI 
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<3. FLOWCHART 



YES 



I 

G.E. = 1.0 




5 


1 


calculate 
1 nduced 
power 


3 


i 


compensate 
for tTp 
1 osses 


1 


f 


compensate 
fo G.E. 







return 
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e. PROGRAMS AND SUBROUTINES USED 



♦ W8S ' 

f. PROGRAM LISTING 



81 ♦LBL "PIT’ 


33 .788 




63 SORT 


0.2 RCL 14 


M * 




64 STO 28 


03 2 






65 Xt2 


84 * 


36 RCL 17 


*' 


66 1/X 


85 SORT 


■5i ■? 




67 RCL 32 


0b RCL 86 


38 YtX 




68 * 


87 / 


39 ~i 1276 




69 2 


88 CHS 


48 * 




78 / 


99 1 


41 + 




71 CHS 


18 * 


42 .5147 

43 + 




72 STO 34 


1! STO S5 




73 Xt2 


12 RCL 89 


44 STO 88 




74 1 


13 RCL 95 


45 (iTG 96 




75 + 


14 2 


46+lBL 85 
47 1 




76 SORT 


15 * 

16 / 




77 RCL 34 

78 * 


17 STO 17 

18 1.5 


48 STO 88 




79 SORT 
88 RCL 18 


19 XOY 


49*LBL 86 




81 * 


28 X>Y? 

21 GTO 85 


58 RCL 25 
51 Xt2 




82 RCL 28 

83 * 


22 RCL 17 


52 STO 32 




84 558 


23 1.3432 

24 * 


53 RCL 18 

54 2 

55 / 




85 / 

86 RCL 15 


25 RCL 17 




87 * 


26 Xt2 


56 RCL 11 

57 / 

58 PI 

59 / 




88 RCL 88 


27 -1.4569 

28 * 




89 * 

99 STO 16 


29 + 




91 *PI=" 


38 RCL 17 


68 RCL 85 




92 ORCL X 


31 3 


b 1 s il 

62 / 




93 AVI EH 


7.? YfX 




94 STOP 

95 END 
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6. Profile pow sr 



a. PURPOSE 

This subroutine calculates the profile power required for 
forward, straight and level flight In terms of horsepower. 

b. ASSUMPTIONS 

Steady flow through the rotor sustem. 

C. EQUATIONS. 

P o - a CdO p A qV ; (1 ♦ 4.25u 2 ) (hp) 

4400 

d. FLOWCHART 



start 



x 

calculate 
prof 1 1e 
povyer 




e. PROGRAMS ANO SUBROUTINE USED 
" WBS " 
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f. PROGRAM LISTING 



01*181. "P01 


02 


RCL 19 


03 


RCL 07 


84 


* 


95 


RCL 11 


96 


* 


87 


RCL 85 


08 


Vf;0 


99 


* 


10 


PI 




¥ 


12 


RCL 13 


13 


3 


14 


YtX 


15 


* 


16 


8 


i 1 

i : 


/ 


18 


STO 21 


19 


RCL 24 


29 


d r - iC - 

? « Lv 


21 


* 


00 
u L 


1 


«57 


+ 


24 


RCL 21 


25 


* 


26 


550 


27 


/ 


28 


STO 21 


29 


°P0=* 


30 


ARCL X 


34 


fiVIEW 


72 


STOP 


TO 


END 
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7. 



Parasite Power 



a. PUPOSE 

This subroutine calculates the parasite power required in 
forward, straight and level flight. 



b. ASSUMPTIONS 

Steady flow through the rotor system. 



c. EQUATIONS. 

PP - 

d. FLOWCHART 






1100 



(hp) 




I 



calculate 

parasite 

power 






store PP 


in 


R28 


! 


t 



return 



S. PROGRAMS ANO SUBROUTINES USED 
« yBS » 
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f. PROGRAM LISTING 



01+LBL “PPl * 
02 SCI 11 
S3 RCL 26 

0.4 * 

05 .5 
86 * 

07 RCL 25 
98 3 
»9 YtK 
18 * 

11 558 

12 / 

13 STG 28 

<4 -pp=* 

15 RRCL X 

16 RVIEw 

17 STOP 

18 END 
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3* Maxim um Forward 7 e lo c i 7 



a. Purpose 

This subroutine calculates the power-1 Iml ted maximum speed 
of the specified helicopter. 

b. ASSUMPTIONS 

a. The power- 1 1ml ted maximum velocity may be estimated 
by neglecting the variatin of Induced power and profile 
power with speed. 

b. Power required to hover is approximately equal to 
power required fo maximum speed. 

c. Steady flow through the rotor system. 

c. EQUATIONS. 

VMAX - v, C-L ] 1/5 

f/A „ 



V 



I 



2 p A 



3 



1/2 
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D. FLOWCHART 
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e. PROGRAMS AND SUBPOUF INE 5 USED 



* W8S * 



f . 



PROGRAM LISTING 



8i*LBL •VHSX' 

82 RCL 26 

83 PI 



84 / 

85 RCL 85 



86 Xt2 



88 i/X 

89 4 
18 * 

11 .3333334 

12 YtX 

13 RCL 28 

14 * 

15 STO 81 

16 1.6884 

17 / 

18 ’YHfiX=* 

19 QRCL X 
28 dVIEH 

21 STOP 

22 RCL 81 

23 RCL 25 

24 JOY? 

25 GTO 12 

26 X<=Y? 

27 GTO 13 



28+LBL 12 
29 *VF > YHflX* 
38 aviE« 

31 YES *CHG* 

32*LBL 13 
33 END 
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9. Best E r.dura ace V = loc i ~ v 

A. PURPOSE 

To calculate the value of velocity corresponding to 
minimum power (l.e., best endurance velocity and/or best 
rate of climb). 

B. ASSUMPTIONS 

1. Steady flow through the rotor system. 

2. The variation of profile power with forward velocity 



Is negl Iglble 



C. EQUATIONS 




D. FLOWCHART 
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e. PROGRAMS AND SUBROUTINES USED 



’ WBS ’ 

f. PROGRAM LISTING 

81*'.Sl = VBE - 

•a; rcl 85 

i?3 sfS 
*4 PI 

85 * 

86 3 

87 / 

88 RCL 26 

89 / 

18 S8RT 

11 RCL 18 

12 * 

13 RCL 11 

14 / 

15 RCL 85 

16 m 



If PT 

iC : * 

19 / 

28 SORT 

21 1.68894 

22 / 

23 -VBEK=" 

24 ORCL X 

25 RVIEM 

26 STOP 

27 EHB 
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10 . 



Srall Or. S3 ~. Velocity 



a . PURPOSE 

This subroutine gives the user of the HP 41-CV an Initial 
approximation for the velocity at which the retreating blade 
angle of attack Is approximately equal to the Is approximately 
equal to the static stall angle of the rotor blade. 

b. ASSUMPTIONS 

1. Steady flow through the rotor system. 

2. Stall onset velocity Is approximately equal to the 



velocity 


for 


best range (I.e., minimum P/v) 


EQUATIONS. 






P 


= 


(-4.25 aCdo A q V t )/(4f) 


r 




(-a Cdo A V 3 )/(4f) 

0 T 


a 


= 


-1/3 * t 


b 


— 


1/27 * (2 * p 2 + 27 * r) 


A 


a 


C -b/2 + (b 2 /2 + a 3 /27) 1/2 ] 1/3 


B 


a 


[ -b/2 - (b 2 /2 + a 3 /27 ) 1/2 ] 1/3 
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d. flowchart 




e. PROGRAMS ANO SUBROUTINES USED 
" WBS " 



33 







f. PROGRAM LISTING 





47 RCL 31 


93 Xt2 


01+LBL "VS" 


48 Xt2 


94 SORT 


92 FS? 03 


49 4 


95 .891 


03 GTD 10 


58 / 


96 XOY 


94 RCL 19 


51 RCL 38 


97 xor? 


05 RCL 97 


52 3 


98 GTO 11 


99 * 


53 YtX 


99 RCL 31 


97 PI 


54 27 


108 RCL 41 


98 * 


55 / 


191 X/Y? 


89 RCL 85 


56 + 


102 GTO 12 


18 m 


57 SORT 


163 RCL 41 


11 * 


58 STO 33 


184 RCL 31 


12 4 


59 RCL 31 


105 X>Y? 


13 / 


69 CHS 


186 GTO 13 


14 RCL 26 


61 2 




15 / 


62 •' 


187+LBL 12 


16 STO 83 


63 + 


108 RCL 25 


17 RCL 13 


£4 . 333334 


109 5 


18 * 


65 YtX 


118 - 


19 4.25 


66 STO 32 


111 STO 25 


29 * 


67 RCL 31 


112 GTO -PGH 


21 CHS 


68 CHS 




22 STO 81 


69 2 


113*LBL 13 


23 RCL 13 


78 / 


114 RCL 25 


24 3 


71 RCL 33 


115 5 


25 YtX 


72 - 


116 + 


26 RCL 83 


73 .333334 


317 STO 25 


27 * 


74 YtX 


118 GTO -OGH 


28 CHS 


75 STO 34 




29 STO 82 


76 RCL 32 


119+LBL 11 


39 27 


77 + 


128 CF 83 


31 * 


78 STO 02 


121 RCL 41 


32 RCL 91 


70 .25 


122 R-B 


33 3 


88 * 


123 *iSTflLL= 


34 YtX 


81 RCL 92 


124 QRCL X 


35 2 


82 + 


125 P.VIEH 


36 * 


83 STO 25 


126 STOP 


37 + 


84 SF 83 


127 RCL 25 


38 27 


85 GTO -9GN* 


128 1.68894 


39 / 




129 / 


48 STO 31 


86*LBL 10 


138 -YS=- 


41 RCL 81 


37 RCL 27 


131 flRCL X 


42 Xt2 


88 RCL 12 


132 RVIEH 


43 -3 


89 / 


133 STOP 


44 i/X 


99 STO 31 


134 END 


45 * 


91 RCL 41 




46 STO 38 


92 - 
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11. I nf low Rat io 
(a) PURPOSE 

This subroutine calculates the ratio of the net velocity up 
through the rotor system to the tip speed. 



(b) ASSUMPTIONS 




Steady flew 


through the rotor system. 



(C) EQUATIONS 




A 


= - / CT / 2 


A 


= - .CT. .... + u Tan a 

2 / a 2 + y 2 3 



D p = ( PP* 5 50) /V F 

a 3 = -Tan (D p /W) 



85 



d 



FLOWCHART 



yes 
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e. PROGRAMS AND SUBROUTINES USED 



* WBS * 



f. PROGRAM LISTING 

6ULEL “LfiHB” 



82 


RCL 


14 


93 


•5 

L. 




84 


/ 




85 


SORT 


86 


CHS 




87 


STO 


63 


98 


,1 




89 


RCL 


22 


18 


X<=Y? 


11 


GTO 


66 


12 


RCL 


28 


13 


558 




14 


* 




15 


RCL 


25 


16 


/ 




17 


RCL 


18 


18 


/ 




19 


STO 


23 



39 


RCL 


83 


48 


- 




41 


Xt2 




42 


SORT 


43 


.88885 


44 


X>Y? 


45 


GTO 


82 


46 


RCL 


88 


47 


STO 


83 


48 


GTO 


81 


49*LBL 


02 


58 


RCL 


88 


51 


STO 


83 


52+LBL 


86 


53 


END 





20*LBL 01 

21 RCL 14 

22 CHS 

23 RCL 63 

24 m 

25 RCL 24 

26 + 

27 SORT 

28 2 
29 * 

38 / 

31 EHTERt 

32 RCL 23 

33 R-B 

34 TSH 

35 RCL 22 

36 * 

37 - 

38 STO 68 
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12 . 



An^ll Cf Attack §5 90 D earae s 



a. PURPOSE 

This subroutine calculates the angle of attack at the 
azimuthal position of 90 degrees. 

b. ASSUMPTIONS 

a. Steady flow through the rotor system. 

b. Blade oscillations are periodic In nature. 

c. Only first harmonics of flapping are necessary for 
calculating angle of attack. 

d. The thrust vector passes through the C.G. 

e. Only uniform twist of the rotor blade Is possible. 



c. EQUATIONS 



a90 



e + e 
0 2 



X 

t ( 1 + y ) 
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d. FLOWCHART 
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e. PROGRAMS AND SUBROUTINES USED 



’ W8S ‘ 
" a 90 ** 



f. PROGRAM LISTING 



81+LBL "i90° 



82 


RCL 33 


83 


RCL 32 


84 


+ 


85 


RCL 29 


66 


+ 


87 


ENTERt 


88 


RCL 83 


89 


1 


!0 


RCL 22 


1! 


+ 


12 


/ 


13 


+ 


14 


STO 42 


15 


R-B 


16 


■498B= 


17 


0RCL X 


18 


OVIEW 


19 


STOP 


28 


END 
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13. 



C ompr es sib ility Pow sr 

a. PURPOSE 

This subroutine calculates the power required due to com_ 
presslblllty on the main rotor sustem In forward .straight 
and level flight In terms of horsepower. 

2. ASSUMPTIONS 

a. Steady flow through the rotor system. 

b. The compresslbll Ity losses can be expressed as a function 
of the amount by which the drag divergence Mach nunber Is 
exceeded at the tip of the advancing blade. 

c. EQUATIONS. 



M 



V f t - .Vi. 



tip 



e 




M 



M 



crlt 



0.06 



t 



M 



crtt 



- 0.71 - 2.3 * a90 



C 



pc 



* CO. 012 * AM ♦ 0.1 * (AM ) S ] 

d d 



p 



. C * p * PI * R* * V s (HP) 

pc t 



550 
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d. FLOWCHART 
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f. PROGRAM LISTING 



81+LBL "CPC" 


37 RCL 19 


82 RCL 13 


38 * 


83 RCL 25 


39 STO 34 


84 + 


48 RCL 11 


85 RCL 43 


41 * 


86 / 


42 RCL 85 


07 STO 3! 


43 X?2 


98 RCL 42 


44 * 


89 RCL 12 


45 PI 


18 * 


46 * 


11 .113 


47 ENTER? 


12 * 


48 RCL 13 


13 CHS 


49 3 


14 RCL 44 


58 Y?X 


15 + 


51 * 


16 STO 32 


52 558 


17 CHS 


53 / 


18 RCL 31 


54 STO 18 


19 + 


55 GTO 10 


28 .86 




21 - 


56+L8L 89 


22 STO 33 


57 8 


23 8 


58 STO 18 


24 XOY 




25 XOY? 


59*LBL 18 


26 GTO 89 


68 ■PC- 


27 RCL 33 


61 fiRCL X 


28 3 


62 PVIEH 


29 YfX 


63 STOP 


38 .1 


64 END 


31 * 




32 ENTER? 




33 RCL 33 




34 .812 




35 * 




76 + 
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14. Stall P owe r 



a. PURPOSE 

Thfs subroutine estimates the additional power required in 
forward, straight and level flight due to retreating blade 
stall. Additionally, this subroutine calculates a stall 
correction factor, k that corrects for the special case. 

s 

of inboard stalling. 

b. ASSUMPTIONS 

a. Steady flow through the rotor system. 

b. The section drag coefficient at stall jumps approxi- 
mately 0.08 at stall onset. 

c. The stalled area is symmetric about the 270 degree 
azimuthal position. 

d. For all airfoils considered, the static stall angle is 
approximately 12.5 degrees. 



c. 



EQUATIONS 

-B 

*1 



2 



* e 

f 



s 1.0 



where, 

B * -u * e - f 

• T 

T - 0.218166 - e ♦ e 

o ? 

8 /2 * e + X 

k =» -( ) 

1 - X 

c - r u >’* <1 - *1 * (1 

D» 



P = C * p * PI * R 2 * V 3 

S D» T 







A/2 



] 
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d. FLOWCHART 




e. PROGRAMS AND SUBROUTINES USED 

* was * 

* SD " 

" CT * 

" W " 

* PP1 " 

" LAMB " 



95 



f. PROGRAM LISTING 







72 * 


81*LBL -CPS’ 


31*LBL 82 


73 558 


02 PS? 0! 


32 *KS=* 


74 / 


93 GTO 86 


33 ORCL X 


75 STO 23 


04 RCL 38 


34 SVIEH 


76 FS? 82 


05 RCL 08 


35 STOP 


77 GTO 87 


86 + 


36 SP 02 




87 2 


37 RCL 15 


78*181 86 


88 / 


38 24 


78 9 


85 1 


35 / 


88 STO 23 


10 XOY 


40 PI 




11 X<=Y’ 


41 / 


81+LBL 87 


12 GTO 81 


42 ENTER! 


32 -PS=" 


13 1 


43 1 


83 flRCL X 


14 3T0 23 


44 RCL 22 


84 flVIEH 


15 GTO 82 


45 - 


85 STOP 




46 X!2 


86 END 


!6*LBL 01 


47 * 




17 RCL 31 


48 1 




18 2 


48 RCL 38 




18 / 


50 - 




29 RCL 25 


51 * 




21 / 


52 1 




22 RCL 38 


53 RCL 38 




23 + 


54 m 




24 CHS 


55 - 




25 ENTER! 


56 SORT 




26 1 


57 * 




27 RCL 38 


58 RCL 23 




23 - 


59 * 




25 / 


60 STO 23 




38 STO 23 


61 RCL 13 






62 3 






63 VtX 






64 RCL 23 






65 * 






66 PI 






67 ♦ 






68 RCL 05 






69 Xt2 






78 * 






71 RCL 11 
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15. Total Power 



a . PURPOSE 

This subroutine calculates the total Dower required In 
forward, straight and level flight, to Include stall and 
compressibility power In terms of hosepower. 

c. ASSUMPTIONS 

Power losses, such as transmission and cooling, can be 
Ignored. 



c. EQUATIONS. 

P =* P + P + P + P + p 

T f o p 3 M 



d. FLOWCHART 




0 



recall P 
power " 


3 




sum a 1 1 
power 
required 


3 


t 



return 
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programs and subroutines used 



e . 



* WBS * 

f. PROGRAM LISTING 



®i*LBL "FT 1 ■ 

82 RCL 28 

83 RCL 16 

84 + 

85 RCL 21 

86 + 

87 RCL 23 

88 + 

89 RCL 18 
18 + 

11 -PT=- 

12 ORCL K 

13 SVIEM 

14 STOP 

15 EHH 
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16 . 



D en s i t v / 5o n i c V s lcci tv 



a. PURPOSE 

To generate values of air density and sonic velocity. 

b. ASSUMPTIONS 

(1) Geopotentlal altitude (H), and geometric altitudes are 
equal below 20.000 feet (actual AH = 29 ft). 

(2) In the troposphere the standard temperature lapse 
rate Is -3.57 *F per 1000 feet. 

c. EQUATIONS 

© * T/T - (1.0 - 6.8753 E-06 * H ) 

SSL 

SVEL = SVEL l * SORT (9) 

RHO = 0.0023769 * (1.0 + HTH * (-.02875 + 0.000275 * HTH) 



where 



HTH = DA/1000 



d. FLOWCHART 
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PROGRAMS AND SUBROUTINES USED 



e . 



• V8S ’ 

f. PROGRAM LISTING 



01+LBL 'BEN" 
92 RCL 45 
93 6,875 E-86 
04 * 

95 CHS 

96 1 

97 + 

98 SORT 

99 1116.89 
10 * 

11 STO 43 

12 RCL 45 

13 1099 

14 / 

15 STO 33 

16 .080275 

17 * 

18 -.92875 

19 + 

29 RCL 33 

21 * 

?V - 

23 + 

24 ,9823769 

25 * 

26 STO I! 

27 END 



ICO 



17 . 



a. PURPOSE 

This subroutine Is used to expedite the changing of up 
to five of the Input parameters whenever a design restraint 
Is exceeded and/or at the end of the main program. 

b. ASSUMPTIONS 

The primary parameters which win require changing are: 



a. 


Forward velocity, VF; 


(Kts) 


b. 


blade twist, TW; 


(rads) 


c. 


1 1ft curve slope, a; 


(per rad) 


d. 


tip velocity, VT; 


(ft/sec) 


e. 


weight, W; 


(lbs) 



c. EQUATIONS 

VF (ft/sec) - VF ( Kts) * 1.66894 
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d. FLOWCHART 




set flags 
04, 2 f 



go to INO 
R31 * 



DISPLAY: 












VF(Kts) 


TW 


a 


VT 


W 



prompt: 

VF(Kts) 



change Kts to 
ft/sec & store 

_ aas. 



prompt: 
TWIST =? 
^tnrg P?q 



prompt: 
vt ~ ? 
store .EL- 



prompt: 
a - ? 
.store -3._ 



prompt: 
W » ? 
ktore.RlQ. 



clear 
flag 27 

1 



NOTE: * returns and 
reruns main 
program wfth 
new data. 
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e. PROGRAMS AND SUBROUTINES USED 



* WBS • 

f. PROGRAM LISTING 



0i*LBL "CNG* 

82 "AGN" 

83 ftSTO 31 

84 CF 84 

85*LBL 86 
86 "CHANGE'* 

87 PROMPT 

88 X=8? 

89 GTO 87 
18 3F 84 
IS SF 27 

22 *VF TH a VI H* 
13 PROMPT 



31+LBL B 

32 "VT=?- 

33 PROMPT 

34 STO 13 

35 GTO 85 



36*LBL E 

37 "U=?" 

38 PROMPT 

39 STO 18 



48*LBL 85 

41 CF 27 

42 GTO 86 



14*LBL A 
15 *VF<KTS>=?" 

16 PROMPT 

17 1.68894 

18 * 

19 STO 25 
28 GTO 85 



43*L8L 87 
44 FS? 84 
45 GTO INS 31 
46 END 



21*LBL 8 
22 "THIST=?" 

23 PROMPT 

24 STO 29 

25 GTO 85 



26+LBL C 

27 *a=?" 

28 PROMPT 

29 STO 68 
38 GTO 85 
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D. HP 41-CV SAMPLE OUTPUT 

Depending cn the value cf forward velocity there a 
three possible types of output that the HP 4 1-cv compuc 
program is capable Of producing. Examples of -.he thr 
possible cases are listed below. 

CASS 1 : Forward velocity, VF, is less than seal 

onset velocity. 



PI=?9i. 751153 
P0=224. 504453 
PF-8. 000688 

¥S= 137. 762848 
V8EK = 76. 069750 
VtfflX=163. 448338 

CYCL1C=8. 888806 
C0LL=17. 896987 

2278=4.388790 

2985=4.888799 

PS=8. 888068 

?M=0. 089860 



P7=l ,816.255686 
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CASS 2 



* h a n 



: Forward velocity, V F, i 

stall onset velocity. 



P!=i 13.87 
P0=352. 43 
P?=784.24 

¥$=137.76 
¥BEK=76. 86 
¥HfiK=!63.45 



CYCLIC--?. 39 
COLL =28. 76 

4278=17.75 

4?8D=-i.83 

M$=8.74 

X0=1.24 

R $= 8.?6 

P$=148.90 

?*=253.71 



PT=1,573.16 



1 05 



greater 



CAS 2 3 



Forward velocity, VF, is greater than, 
the maximum forward velocity possible. 



PI=107. 18 
P0=368.92 
PP=8«4.72 

V$=137,76 
VBEK=7M6 
VtfRX-163. 45 
VF > VHflX 

CYCL I C=- 1 8 . 55 
COLL =22. 88 

4270=19,87 

i98B=-1.23 

X$=S,73 

X0=1.53 

KS=1.88 

PS-158,18 

?*=289.23 



PT= 1.-786. 15 
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A. IBM 3033 PROGRAM DOCUMENTATION 

This program calculates the power required fly a 

helicopter in forward, straight and leveL high speed flight. 
In the sections mat follow all inputs and implementation 
requirements are specified. The problem solving methodology 
used is as described in Chapter II. 
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3. INPUT DATA REQUIRED 

In calculating the performance of a helicopter it is 
necessary to define a set cf force conditions, environmental 
conditions, and physical conditions. The force conditions 
required are the aircraft gross weight (W) in pounds, the 
maximum lift coefficient (CLMAX) , and ths the coefficient of 
list at zero angle of attack (CDO) . The environmental 
cor.ditons include, forward velocity (VFK| in knots, speed of 
sound (SVEL) in feet/sec, rotor height above the ground (H) 
in feet, and air density (RHO) in Ib-ft /sec . For sake of 
user simplicity, the the speed of scund and air density are 
generated within the program and are comparable to the 
values found in a standard atmosphere table. Finally, the 
physical conditions required are the rotor radius (R) in 
feet, tip velocity (VT) in ft/sec, numoer of blades (b) , 
main rotor chord (C) in feet, flat plate area (FPA) in ft , 
geometric twist of the rotor (TWIST) in radians, and airfoil 
lisft curve slope (CLA) in per radians. 
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C. HELICOPTER SAMPLE DATA 

Table C.1 below illustrates the format required when 
inputing data. 





TABLE C.l 


HELICOPTER SAMPLE IN POT 


DATA 
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CDO 
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£A 
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0.0 
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1.4 


2. 


7. 




1000 
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3M PROGRAM FLOWCHART 
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profile 
power 
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S. IBM PHCG RAM LISTING 

This section contains the listing of the IBM 
computer program developed in - h is report. 
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A AAAAA A A A A A 



c 



ENO 



MAIN ROTOR HEIGHT ABOVE GROUND. .. ',012.o/, 

AIR DENSITY IRHGI *,F12.o/, 

SONIC VELOCITY *,012.6/, 

BLACc GEOMETRIC TWIST............ *,(ji2»6//» 

TYPE A IRFCI L • ,012.6/ , 

MAXIMUM 2-0 LIFT COEFFICIENT..... ',012.6/, 
2- D STATIC STALL ANGLE (AMAX).... ',012.7/, 

LIFT CURVE SLOPE (/RAD) ‘,012.6/, 

2ERC-LIFT DRAG CO E EPIC I ENT F 12. 6 / , 

CRITICAL MACH NO ( FCR CL = 0) . . • , Fi2. 6 / / 
INITIAL FORWARD VELOCITY (XT).... ',012.6/) 
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SAMPLE OF IBM COMPUTES OUTPUT 



This 


section contains 


an example run of the 


IBM computer 


orogran , 


utilizing AH1-J 


Cobra data, scarring 


at a forward 


velocity 


cf 120 knots and 


terminating a: VMAX. 
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HIGH SPEED FORWARD FLIGHT ANALYSIS 



INPUT DATA 



RADIUS 

MAIN ROTOR CHCRC 

NUMBER OF MAIN FQTOR BLADES 

AICRAFT GROSS WEIGHT..' 

ROTOR TIP VELOCITY 

HORIZONTAL FLAT PLATE AREA 

NUMBER OF ENGINES IN HELICOPTER.. 

DENSITY ALTITUDE 

MAIN ROTOR HEIGHT ABOVE GROUND... 

AIR DENSITY (RHCJ 

SONIC VELOCITY 

BLADE GEOMETRIC TWIST 

TYPE AIRFOIL 

MAXIMUM 2-D LIFT COEFFICIENT 

l 2 T?t s mi SlfiW 

ZERO-LIFT CRAG COEEFICIENT 

CRITICAL MACH NC (FOR CL = 0).. 

INITIAL FORWARD VELOCITY (KT).... 

COEFFICIENT CF THRUST.... 

INDUCED VELOCITY 

DISC AREA 

SOLIDI TY 

TIP-LOSS FACTOR 

GROUND EFFECT RATIO 



22.0000 

2.25000 

2.00000 

10612.0 

738.000 

17.0000 

2.00000 



1000.00 
1000.00 
0. 00230 
1113.04 
-.174500 

7.00000 

1.40000 

■mm 

0. 01075 
0. 72000 



9 



0 

0 



12 0. COO 

• 554153 E -02 
33. 8469 

15 22 .62 

• 65 0197 E-01 
.947362 
1.00000 



£***#$*#*****# *** #**#**#****#**#*****$#***:»***:>#* 
FORWARD VELOCITY IN KNOTS = 120. OCO 



PARASITE DRAG = 806.256 

DISK PLANE ANGLE CF ATTACK = -.759760E-01 



DYNAMIC PARAMETERS 



A 

A 

A 

A 

T 

T 

T 

T 

T 



11 

12 

13 

14 

1. 



4 

NFLC'W RATIC 



LAM 6 CAi 



0.625400 
0. 806926 
0. 573338 
1.175435 
0. 4676 02 
0.319142 
0.216296 
0.125826 
— 0 . 0403 8 1 



STALL PGWER CALCULATIONS 



R T 1C6773E-0 i 

IN80ARD STALL CORRECTION FACTOR 0 

STALL POWER COEFFICIENT 



ANGLE OF ATTACK CALCULATIONS 



LONGITUDINAL CYCLIC ANGLE -0. 109613 

LONGITUDINAL COLLECTIVE ANGLE... 0.314954 

A LPHA ( 90 ) (DEG ) -0.043036 

A LPHA (270) (DEG) 12.511999 



HIGH SPEEC MACH EFFECTS 

ADVANCING ELADE TIP MACH NUM3ER... .345135 

CRITICAL MACH NUMBER 720543 

DRAG DIVERGENCE MACH NUMBER 645919E-01 

COMPRESSIBILITY POWER COEFF 521491E-Q4 



444444444 ****** 4444444444444444444444444444444444**444 



PCWER REQUIRED 

44 44444444444444444444444 *44444444*4*4444444444*444444 



INDUCED PCWER = 
PROFILE PCWER = 
PARASITE POWER = 

COMP RES I B IL IT Y POWER 
STALL power = 

TOTAL POWER REQUIRED 



1 51. 545 
296.464 
29 7. 103 
134. 001 

.0 

379.112 



4444444444 4444 444 4 4*4 4*4 44 4# 4 414c# 4c 4:4# 44 4**4 4#*# 44= * 



FCRWARO VELOCITY IN KNOTS = i30.000 

4444*4444444444444 444 4 44 44444 444* 44444 444 4 44 44 4*4 



PARASITE DRAG = 946. 234 

DISK PLANE ANGLE OF ATTACK = -.6916o4E-01 



DYNAMIC PARAMETERS 



All..... 

A 12 

A 13 

A14.. 

T 1 

T 2 

T 3 

T 4 

INFLOW RATIO ( LAMBDA) 



0. 680172 
0.880377 
0. 625882 
1.207472 
0.470875 
0. 325344 
0.221234 
0. 136796 
-0. 044668 



STALL FGnER CALCULATIONS 

RT...... •••••••••• . . . • 75341 cc— 02 

INBOARD STALL CORRECTION FACTOR 0 

STALL POWEF COEFFICIENT 0 



ANGLE OF ATTACK CALCULATIONS 



LONGITUDINAL CYCLIC ANGLE -0.122350 

LONGITUDINAL COLLECTIVE ANGLE... 0.326139 

ALPHA ( 90 I (QEG) -0.291443 

A LPHA ( 27 0 ) (DEGi 13.728103 



HIGH SPEEC MACH EFFECTS 

ADVANCING eLAGE TIP MACH NUMBER... 

CRITICAL MACH KUMBER 

DRAG DIVERGENCE MACH NUMBER 

COMPRESS I El c IT V POWER COEFF....... 



. 860309 
.723294 
•770157 E- 01 
.63Q606E— 04 



494999949444494499949944999999994949999999999999999999 

POWER RE GUI RED 

4 4 4444 4444 4 444 4 444 4 444 944444444444 99999999949499999994 



INDUCED POWER = 
PROFILE POWER * 
PARASITE PCWEP = 
COMPRESI 3IUT Y POWER 
STALL POWER = 

TOTAL POWER REQUIRED 



139.900 
308. 957 
377.740 
162. 039 

.0 



938 • 63o 



4 9999 994 4 44499 494 9 4 49 449 944494449 44449 499 4 4999 499 

FORWARD VELOCITY IN KNOTS = 140.000 

4 4444994 994449 4999 499 449 94 9994999 4999 9 499 9 4949 499 



PARASITE DRAG = 1097.41 

DISK PLANE ANGLE OF ATTACK = -.103412 



DYNAMIC PARAMETERS 



All 0.735630 

A 12 0. 55656 3 

A 13 0. 679055 

A 14 . . •••.. .... 1 . 2426 2 5 

T 1 0.474411 

T 2 0.332042 

T 3 0.224407 

T 4 0. 147338 

INFLOW RATIC { LAMBCA) -0.050098 
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STALL POWER CALCULATIONS 



RT . -.389701E-02 

INBOARD STALL CORRECTION FACTOR 0 

STALL POWER COEFFICIENT 0 



ANGLE OF ATTACK CALCULATIONS 

LONGITUDINAL CYCLIC ANGLE -0.13657 6 

LONGITUDINAL COLLECTIVE ANGLE... 0.335933 

A LPHA ( 90 ) (DEG) -0.520492 

A LPHA (270) (DEG) 15.129103 



HIGH SPEEC MACH EFFECTS 

ADVANCING ElADE TIP MACH NUMBER... 

CRITICAL MACH NUMBER 

DRAG DIVERGENCE MACH NUMBER 

COMPRESSIBILITY POWER COEFF 



. 875464 
.725882 
. 89601 2 E-Oi 
. 7453 7 2E - 04 



444444444444444 444=4 444 44444444444444444444 4 44 4 4 4jp:?444 4 

POWER REQUIRED 

4444444444444444444 444 444 4 4444 444 4 444 4 4444 4 44 4 4 444 444 * 

INDUCED POWER = 129.913 

PROFILE POWER = 322.449 

PARASITE FCwEf = 471.789 

COMP RES I 8 I LI T Y POWER = 191.657 

STALL POWER = .0 

TOTAL POWER REQUIRED = 1115.81 



4 4444 4 44 4 4 4 4 4 4 4444 444 444 4 4 44 4 4444 444 4 4444 4 44 4 4 444 



FCRWARC VELOCITY IN KNOTS = 150. OCO 

4 4444 444 4 44444 444 4 444 444 44444444 4 44444 444 4 44 44 444 



PARASITE DRAG = 1259.78 

DISK PLANE ANGLE OF ATTACK = -.118713 



DYNAMIC PARAMETERS 



All 

A 12 

A 13.. 

A 14.. 

T 1 

T2. 

T 3 

T 4 

INFLOW RATIC ( LAMBDA) 



0. 791349 
1. 03*1 7C 

0. 724801 

1. 2810 51 
0. 478208 
0. 339237 
0. 227315 
0.159103 

-0 • 0567 2 0 
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STALL P GWER CALCULATIONS 



R T 



.2307 77E-02 



VALUES FOP DETERMINING STALL 



GAMMA 

CS 

BS 

XS 

XO 



-C. 2 6 459 3 SO 
-C. 14755023 
0.32449627 

0. 7921403 

1 . C674362 



CQEFFICI ENT 



INBOARD STALL CORRECTION FACTOR 6622 18C94 

STALL POWER COEFFICIENT .31202595E-0* 



ANGLE OF ATTACK CALCULATIONS 

LONGITUDINAL CYCLIC ANGLE 

LONGITUDINAL CCLLECTIVE ANGLE... 

ALPHA (90) ( CEG ) 

ALPHA ( 270 ) (DEG) 



-0. 152545 
0.356377 
-0.73866 7 
16.740753 



HIGH SPEEC MACH EFFECTS 

ADVANCING BLADE TIP MACH NUMBER... 

CRITICAL MACH NUMBER 

DRAG DIVERGENCE MACH NUMBER.. 

COMPRESS I BILIT Y POWER COEFF 



.890658 
. 728348 
.102309 
. 86733 3c- 04 



******************* *** ************ ******************** 
PCWER RE GUI RED 

******************* *** ******************************** 



INDUCED PCwER = 
PROFILE PCWER = 
PARASITE POWER = 
COMPRESI BILITY POWER 
STALL POWER = 

TOTAL POWER REQUIRED 



121.281 
336. 9 41 
580. 279 
223. 009 
80. 1 774 

1341 .69 



************** ******* ************ **************** 
FORWARD VELOCITY IN KNOTS = 160.000 

* ******* ****** ******* **************************** 



PARASITE DRAG = 1433.35 

DISK PLANE ANGLE OF ATTACK = -. 



13 50o3 



DYNAMIC PARAMETERS 



A li 

A 12 

A 13 

A 14 

Tl... . 

12 

T3 

T 4 

INFLOW RATIO ( LAMBDA) 



0. 8489 C 7 
1, 113397 
0. 791449 
1.322899 
0.482267 
0.346927 
0. 221457 
0. 170453 
-0. 064591 



STALL POWER CALCULATIONS 
RT 



.1158525-01 



VALUES 

GAMMA = 

CS 

BS 

XS 

XO 



FOP DETERMINING STALL COEFFICIENT 

-C. 30194968 
-C. 17515445 
0.3 6584556 

0. 739 8 598 

1. 2566751 



INBOARD STALL CORRECTION FACTOR l.OCOOOCOO 

STALL POWER COEFFICIENT 6054921 CE-04 



ANGLE 0 F ATTACK CALCULATIONS 

LONGITUDINAL CYCLIC ANGLE 

LONGITUDINAL COLLECTIVE ANGLE... 

ALPHA (90) (DEG) 

ALPHA(270) (CEO... 



-0. 17C535 
0.375743 
-0. 94936 i 
18.59143 i 



HIGH SPEEC MACH EFFECTS 

ADVANCING ELADE TIP MACH NUMBER... .905332 

CRITICAL MACh NUMBER 73C729 

DRAG DIVERGENCE MACH NUMBER 115102 

COMPRESSIBILITY POWER COEFF 997223E-04 



***#*$.#:******:*:* *** ? *** *** ***.*:$## ** c# ** #$***«** * 

POWER RE GUI RED 



INDUCED POWER = 113.722 

PROFILE PCWER = 352.433 

PARASITE POWER = 7C4. 243 

COMP RESIBILITY POWER = 25o.243 

STALL POWER = 155. 53o 



TOTAL POwER REQUIRED 



1582 .23 



FORWARD V 2LQC ITY IN KNOTS = i63.359 

* *** **«*$:* ***#*4:# *#* **##***:* #**»#**;£# 



PARASITE DRAG = 1494.16 

OISK PLANE ANGLE uf ATTACK = -.140799 



DYNAMIC PARAMETERS 

Ail 

A 12 

A 13 

Alt., . 

T 1 

T 2. • • • 

T3 

T 4. • 

INFLOW RATIC (LAMBDA).. 



stall power calculations 

RT 1559876-01 



VALUES FOF DETERMINING STALL COEFFICIENT 

GAMMA = -C.3 1570 089 

CS = -C. 18555 188 

3 S = 0.38093823 

XS = 0. 7336488 

XO = 1. <493771 



INBOARD STALL CORRECTION FACTOR 1. 00000000 

STALL POWEF COEFFICIENT .61 1 C876CE-04 



ANGLE OF ATTACK CALCULATIONS 

LONGITUDINAL CYCLIC ANGLE....... -0.177082 

LONGITUDINAL COLLECTIVE ANGLE... 0.382947 

ALPHA i 90 ) ( DEG ) -1.019017 

A LPHA (270 ) (DEG)... 19.272018 



HIGH SPEEC MACH EFFECTS 

ADVANCING BLADE TIP MACH NUMBER... .910929 

CRITICAL MACh KUMBER .731516 

DRAG DIVERGENCE MACH NUMBER 119^12 

CGMPRcSSI BILIT Y POwER COEFF. 104241E-03 
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0.868277 
1. 141191 
0. 810842 
1.337756 
0. 483689 
0 • 349o22 
0.232734 
0.1742 9 7 
0.067526 






PCwER REQUIRED 

a*#**## :*:{<* £41 aj$ aja* ^4t*t a* * ap# afr+v###*#*###**^** afa$a 4?*$* * 



INOUCED PChER = 
PROFILE PCWER * 
PARASITE FCWEF = 
COMPRESIBILITY POWER 
STALL POWER = 

TOTAL POWER REQUIRED 



111.338 
357. 861 
749.536 
267. 855 
157.023 



1643 .66 



INITIAL STALL CNSET VELOCITY APPROXIMATION (KTSi •= 76 . 



DIFFERENCE BETWEEN A270 AND A, VAX * .29987IE-03 

ACCEP ATA6L E ERROR = .1000C0E-02 



44444444444444 *##* 4 *# 4 ajaaja aja a*** * 44 ** * 4*4 4 *44 4 44 4 4 4 



ANCLE AT STALL ONSET= 13.9811 

STALL CNSET VELOCITY = 131.913 

VELOCITY MAX E I'DURANC E = 76.0347 

MAXIMUM FCPwARC VELOCITY = 163.359 

4 a? a{aa£aja 4 ^a^c^c aja a} 4 44 4 4 4 4 444 444 4444 4 4444 444 44 444 4 44 44 4 
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G. COMPARISON 0? PROGRAM OUTPUT VS TEST D AT A 

This section compares the output cf the HP41-CV and 13 
3033 computer programs tc actual flight rest data gathers 
at the Naval Weapons 'System Center. 
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POWER REQUIRED VS VELOCITY 
FOR THE AH1-J USING THE HP 41-CV 




Figure c.1 



Power Curves Generated 07 HP41-C7 Program 



POWER REQUIRED VS VELOCITY 
FOR THE AH1-J USING THE IBM 3033 



I 




Figure C. 2. Power Curves Generated by the I3H 3033 Prograa. 
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POWER REQUIRED VS VELOCITY 




Figure G.3» Comparison Between Computer Data and Actual 
Test Plight Data. 
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